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This article investigates the signatures of various models of dark energy on weak gravitational lens- 
ing, including the complementarity of the linear and non-linear regimes. It investigates quintessence 
models and their extension to scalar-tensor gravity. The various effects induced by this simplest 
extension of general relativity are discussed. It is shown that, given the constraints in the Solar 
System, models such as a quadratic nonminimal coupling do not leave any signatures that can be 
detected while other models, such as a runaway dilaton, which include attraction toward general 
relativity can let an imprint of about 10%. 

PACS numbers: 98.80.Cq, 98.80.Es,04.80.Cc 



I. INTRODUCTION 

The growing evidences for the acceleration of the uni- 
verse have led to the formulation of numerous scenarios 
of dark energy (see Ref. El for reviews). In this con- 
struction various routes have been investigated 0] ■ They 
reduce mainly to the introduction of new degrees of free- 
dom in the cosmological scenario, either as matter fields 
with negative pressure or in the gravitational sector. 

The simplest models accounting for this acceleration 
rely on the introduction of a slow rolling scalar field 
whose expectation value varies during the history of the 
universe. These quintessence models |3j are characterized 
by their potential and numerous choices have been dis- 
cussed in the literature EI- Generically this scalar field 
may also couple to matter fields Q. As a consequence, 
the values of the fundamental constants will depend on 
the value of this field and may vary [5j and it will induce 
a fifth force that will be long range if the scalar field is 
light, which is usually the case if it is the source of the 
acceleration of the universe. 

In the simplest extension of the quintessence mod- 
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els, the quintessence field universally couples to the 
other fields. It was realized that the properties of 
this quintessence field were conserved in that situa- 
tion 0, 0, • This implies that we are dealing with a 
scalar-tensor theory of gravity in which the spin-0 partner 
of the graviton is also the quintessence field. From a phe- 
nomenological point of view, extended quintessence mod- 
els are the simplest well defined theories in which there 
is a modification of gravity. Various observational signa- 
tures on the background evolution,the cosmic microwave 
background (CMB) anisotropies [1 U E3 EH El, the 
big-bang nucleosynthesis have been worked out [13( and 
the reconstruction problem was discussed in details |14| . 
In particular, it was shown that the attraction mecha- 
nism toward general relativity E3 still hold Eil in these 
extended quintessence models, a crucial point since they 
have to satisfy sharp constraints in the Solar System 1 171 . 
It was also pointed out that a runaway dilaton [Tsl llfj| 
that does not couple universally, is an appealing models 
with specific signature such as a variation of some fun- 
damental constants and a violation of the universality of 
free fall. 

Scalar-tensor theories are the most natural alternative 
to general relativity 0, |2(J , preserving the universality 
of free fall and constancy of all non gravitational con- 
stants. Gravity is mediated not only by a massless spin-2 
graviton, corresponding to the spacetime metric, but also 
by a spin-0 scalar field. Many theoretical motivations to 
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consider such a scalar partner to the graviton have been 
put forward, particularly in higher-dimensional theories. 
In string theory, the supermultiplet of the graviton con- 
tains a dilaton and moduli fields pi). 

Scalar-tensor theories are well constrained in the Solar 
system and one can try to extend these constraints to 
astrophysical and cosmological scales. In the case of ex- 
tended quintessence, the effect of the scalar field is impor- 
tant particularly in the recent universe when it starts to 
dominate. It was pointed out that the modification 
of the equation of state in the recent universe has strik- 
ing effects on the growth of density perturbation and on 
weak lensing observables such as the convergence power 
spectrum. 

Weak gravitational lensing has now proven to be a 
powerful tool to study large scale structures jU 0, 113 
and to gather information on the nature of the dark 
energy (2(| I^J. In particular, cross-correlation tech- 
niques seem a very promising way to achieve this task 
(see e.g. p^)- Weak lensing can be detected by the de- 
formation of the shape of background galaxies. It was 
recently observed by various groups 29\ and can probe 
the large scale structures of the universe both in the lin- 
ear and non-linear regimes. Recently, the weak lensing 
has been studied in the context of generalized cosmolo- 
gies [3(J, specializing to quadratic non- minimally cou- 
pled models and focusing on the effects on the CMB 
anisotropics, hence just looking after the linear regime 
of structure formation. However, as was shown on vari- 
ous examples p^, [2(| E^J , the most stringent constraints 
on cosmological models, able to distinguish dark energy 
models, arise from the comparison of the linear and non- 
linear regimes. 

In this article, our main goal is to study in details the 
lensing observable, focusing on 2-point functions to start 
with, in the case of general relativity and scalar-tensor 
theories of gravity. In particular, most of the previous 
studies (e.g. Refs. P2, p(| [i^]]) treats the matter power 
spectrum independently of the CMB anisotropics. We 
choose to normalize all our spectra on the CMB at low 
multipole and deduce the lensing observables, both in the 
linear and non-linear regimes, with the same normaliza- 
tion. Independently of any considerations on the law of 
gravity, this tool will be of first importance to deal with 
combined analysis of CMB and lensing data [32 ■ 

Our analysis will be applied to both quintessence and 
extended quintessence models. We will show that the 
modification of the equation of state of the universe leaves 
a detectable imprint on weak lensing observables. Con- 
cerning scalar-tensor theories we will show that, given 
the constraints in the Solar System, many models will 
let very little signatures. This is the case for instance of 
a quadratic coupling. Interestingly other models such as 
a runaway dilaton can have a 10% effect at that is likely 
to be constrained and/or detected (see Fig. IT2)l . Note 
that both conclusions are of interest since it will tell us 
for which class of modifications we have to bother about 
this extension of the law of gravity. 



The article is organized as follows. We first introduce 
scalar-tensor theories of gravity in § [Q] In particular, 
we describe the Einstein and Jordan frames and argue 
that the latter is the one in which observations take their 
standard interpretations. We also recall (§ III Bill C)l the 
standard constraints on these models and discuss briefly 
the properties of gravitational lensing (§ III Dfl . In § IIIII 
we derive the distortion of a geodesic light bundle due 
to large scale structures in a way that is valid for any 
metric theory of gravity, in particular it does not assume 
that gravity is described by general relativity and it is 
valid for scalar-tensor theories. This allows us to define 
in § IIVI the weak lensing observables such as the shear 
power spectrum and the 2-point statistics of the shear 
field. We then describe our numerical implementation 
(§ IIV D|) in details, as well as the mapping to the non- 
linear regime. As a first check, we investigate in § Ivl 
weak lensing in general relativity both for a ACDM and 
quintessence models. Sec. I VII discusses the various ef- 
fects and differences that arise in scalar-tensor theories 
of gravity and we then investigate in § IV11I two families 
of models: a non-minimally coupled scalar field and run- 
away dilaton-like models that include attraction toward 
general relativity. 



II. GENERAL RESULTS ON SCALAR-TENSOR 
THEORIES 

In this article, we focus on scalar-tensor theories of 
gravity described by the action 

1D7TG* J 

+S m [g llv ; matter fields] (1) 

where G* is the bare gravitational constant from which 
we define k„ = 87rG». The coupling function F mul- 
tiplying the Ricci scalar R is dimensionless and needs 
to be positive to ensure that the graviton carries positive 
energy. S m is the action of the matter fields that are cou- 
pled minimally to the metric, hence ensuring the univer- 
sality of free fall. The metric has signature (— , +, +, +) 
and we work in units in which c = 1. 

A. Field equations 

The variation of action leads to the field equations 

1 2 

+df,ipd„ip - -g^ (d a f) - g^Uitp) 
+V„d 1/ F(<p) - g^DFfa) (2) 
Op = U, v --F tV R (3) 

v^r" = o. (4) 
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Here is the covariant derivative associated to g^ lL , and 
the subscript ", 92" denotes the functional derivative with 
respect to ip. The stress-energy tensor is defined as 

2 5S m 
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Action Q has been written in the so-called Jordan frame 
in which matter is universally coupled to the metric. The 
Jordan metric defines the length and time as measured 
by laboratory apparatus. In the following, we will be in- 
terested in particular in the shape of background galaxies 
and all observations will have their standard interpreta- 
tion in this frame |2(| . 

It is however useful to define an Einstein frame action 
through a conformal transformation of the metric 



9L = F(<p)9nv- 



(5) 



In the following all quantities labelled by a star (*) will 
refer to Einstein frame. Defining the field ip* and the two 
functions A(ip*) and V(ip*) (see e.g. Q) by 



dip* 
d<p 

A(<p 



3 
4 

= F 



d\n F(tp) 



dip 
1/2 M 

2V(ip*) = U(p)F- 2 (p) 
the action (JJJ reads as 

^-Jd^V^lR* 

+S m [A 2 (ip*)g* ;ip] 



2F(<p) 



S 



^9* L 'dp<p*d u tp* 



(6) 

(7) 
(8) 



(9) 



The kinetic terms have been diagonalized so that the 
spin-2 and spin-0 degrees of freedom of the theory are 
perturbations of g* pv and ip* respectively. In this frame, 
the field equations take the form 



G* 



8ttG*T* 



a*ip* = V, v . - A-kG^*)T; v 9^ 



28^*8^* 

v M rr = a(ip*)T; p g^8„ip* 

where we have defined the Einstein frame stress-energy 
tensor 

SSrn 



9% - 2g;„V (10) 

(11) 
(12) 



r ^9~*^9lv 

related to the Jordan frame stress-energy tensor by T*„ 



A 2 T )1U . The function 



a(ip*) = 



d\nA 



dip* 



(13) 



characterizes the coupling of the scalar field to matter 
(we recover general relativity with a scalar field when it 
vanishes). For completeness, we also define 



/%*) 



da 

dip* 



(14) 



TABLE I: Notations in Jordan and Einstein frames. 





Jordan frame 


Einstein frame 


coordinates 


OM) 


{x„U) 


metric 






scalar field 


<P 


<p* 


potential 


U 


V 


coupling 


F 


A 



Note that in Einstein frame the Einstein equations {HS} 
are the same as the ones obtained in general relativity 
with a minimally coupled scalar field. Table [I] summarize 
the notations used in Jordan and Einstein frames. 

From these definitions, we can define an effective grav- 
itational constant as 



GW = 



G* 
F(<p) 



= G*A 2 (ip*) 



(15) 



but, as we shall see, this does not correspond to the 
Newton constant that will be effectively measured in a 
Cavendish- like experiment. 



B. Local constraints 

Deviations of the theory of gravity from general rela- 
tivity are sharply constrained in the Solar System |17| as 
well as by binary pulsars ;33j ■ The constraints are usually 
set on the post-Newtonian parameters |l7j. In t he p ar- 
ticular case of scalar-tensor theories, they reduce fl4l l20j | 
to 



„PPN 



F 2 



F + 2F 2 



<p 



= -2 



a 



PPN 



1 = -- 



1 



or 

„2 



do 



2(1 
1 FF 



2 \2 



dip* 



PPN 



4 2F + 3F 2 dip 



(16) 

(17) 
(18) 
(19) 



Solar System experiments set sharp constraints on the 
values of the PPN parameters (jq pn , /3 ppn ) today. The 
observed value of the perihelion shift of Mercury implies 
the bound |34| 

|27o PPN -/9o PN -l| <3x 10- 3 . (20) 

The Lunar laser ranging experiment |35| imposes 



4/? PPN - 7o PN - 3 = -(0.7 ± 1) x 10- 3 (21) 



and the measurements of the light deflection by Very 
Long Baseline Interferometry |3j| improves the con- 
straint on 7 PPN to 



|7o PPN 



1 = 4 x 10" 



(22) 
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The recent analysis of the frequency shift of radio waves 
to and from the Cassini spacecraft have set the even more 
stringent bound 



7 PPN 



- 1 = (2.1 ±2.3) x 10 



-5 



(23) 



The previous bounds to those obtained in Refs. 36? ] 
were (see Ref. E3 for a review) 



|7o PPN -l| < 2xl(T 3 , 
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PPN 



1 < 6x 10" 



(24) 



These constraints can be transformed to constraints on 
the set (ao,/3o) (see Ref. [33| for a summary). In par- 
ticular, we can note that binary pulsars imply that 
0o > -4.5. The bound (|23J) implies that 



'2 



F 



< 4 x 10" 



(25) 



The local value of the gravitational constant is deduced 
from a Cavendish-like experiment, i.e. from the measure- 
ment of the Newton force between to masses m i and m2, 
as G cav = Fr 2 / m i m 2- H can be shown [2(j that its 
theoretical expression is 



Gr-av — 




(26) 
(27) 



Note that today Eq. 125fl implies that G cav and G e ff do 
not differ by more than some 10 -3 percent, which is any- 
how larger than the accuracy of the measurement of the 
gravitational constant. At higher redshift this difference 
may be larger. 

Current constraints [38[ (see also Ref. Q for a review) 
imply that 



G r 



< 6 x lO^yr- 1 - 



(28) 



Using Eqs. I|13|) and (|14[) . this implies that 

/3o«o 



ao + 



(i + c4) 



|vs,o| < 3 x 10- 12 yr _1 (29) 



All these constraints are local and we will consider in 
this article only models that satisfy them. 



C. Cosmological constraints 

From a cosmological point of view, there are very few 
constraints on scalar-tensor theories. 

CMB observations may in principle give some con- 
straints but they are often degenerate with other param- 
eters, such as ej?. the cosmological parameters (see e.g. 
Ref. d H E3 El for some studies on the CMB im- 
prints of scalar tensor theories). 



More stringent constraints arise from Big Bang nucle- 
osynthesis (BBN). BBN results, and in particular the 
helium-4 abundance, are very sensitive to the weak in- 
teraction freeze-out. This temperature depends on the 
strength of the gravitational interaction which dictates 
the expansion rates. BBN mainly requires that (i) the 
universe is dominated by radiation at the time of nucle- 
osynthesis and (ii) that the number of degrees of freedom 
of relativistic particles does not vary by more than 20% 
with respect to its expected value g T = 10.75. The Fried- 
mann equation (| A6|) for a flat universe reduces to 



U =n Q -a 



1+ % )T 4 



(30) 



with Sg x / 'g r — F /F(z nuc ) — 1, if one neglects the contri- 
bution of the scalar field to the energy density. Requiring 
that Sg r /g r be smaller than 0.2 implies EH that 



0.8 < 



F 




A 2 

nuc 


TP 

± nuc 




A 2 



< 1.2. 



(31) 



Let us emphasize that large value of \Fq /F nuc \ were shown 
to be consistent with the BBN constraints if (3 is 
large enough so that the naive limit l|31(l can be much 
more stringent than a detailed study may show. 



D. A remark on gravitational lensing 

Gravitational bending of light by a single mass M in 
scalar-tensor theory can be derived easily in Einstein 
frame. Since photons are coupled to the gravitational 
field only, they are insensitive to the scalar field. One 
deduces that they are deflected by an angle A8 exactly 
in the same way as in general relativity, that is 



A9 = AG^A 2 (^)M/r , 



(32) 



where A 2 M is the deflecting mass in Einstein frame. G* 
cannot be measured directly and what we know is only 
the value of the gravitational constant l|26|) determined 
by a Cavendish-like experiment today. Masses attract 
themselves due to the exchange of both a graviton and a 
scalar. It follows that 



G cav = G*Ao(± + 2 



G*Aq 



l + 7o 

so that the deflection angle turns out to be 
A0 = 2(l + 7 o)G cav M/r o . 



(33) 



(34) 



In conclusion, the deflected angle is different from the 
prediction of general relativity, not because photons are 
deflected differently but because there is a scalar inter- 
action between the test masses which contributes to the 
value of the effective Newton constant [2^ • 

This argument can be straightforwardly generalized to 
the cosmological context. If the lens seats at a redshift z 
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and if we work in the thin lens approximation then the 
deflection angle is given by 



A9 = 4G»A 2 ( 

-Ions 



(35) 



again, G* cannot be measured and we have access to 



G c 



i G*Aq 
l+7o 



(36) 



where the subscript means that it is determined today. 
It follows that 



A0 = 2(l+ 7o )G, 



cav.O " 



A/lcns A 2 (z) 

r A*(0)' 



(37) 



The mass of the lensing galaxy, Mi cns , cannot be mea- 
sured directly. 

We mainly expect to have two major effects on lensing 
observations: (i) an effect of the background dynamics 
through the angular distances, (ii) a gravitational effect 
due to the z dependence deflection angle. Eventually, 
observations of lenses located at different redshift should 
enable to give an information on A(z)/A(0) as a function 
of the redshift. 

Before we go into a detailed computation, we can 
give an upper bound on the amplitude of this effect. 
From the constraint (|3f |) . we deduce that the factor- 
depending on A in Eq. (pITjl has to lie in the range 
F /F(z nuc ) ~ [0.8,1.2] so that effects larger than 20% 
are unlikely to happen, independently of any model (see 
however Ref. Let us stress that, as we shall see, the 

dominant effect will arise from the growth of the density 
field and not from the two effects we have just mentioned. 

To finish this general discussion let us emphasize that 
we could have added a coupling of the scalar field to the 
electromagnetic tensor of the form 



(38) 



in the action Such a term will induce a variation 
of the fine structure constant (see e.g. Ref. |5j for a re- 
view), that has to be smaller than 10~ 5 between z ~ 3 
and today, and a violation of the universality of free fall. 
Interestingly, such a term will not affect the equation of 
propagation of photons in the eikonal approximation at 
first order and thus weak lensing observables that are 
considered in this article. 



III. LIGHT PROPAGATION IN A PERTURBED 
FLRW SPACETIME 

This section is devoted to the general theory of gravita- 
tional lensing. Our construction is based on the geodesic 
deviation of a bundle of null geodesies and will then be 
valid as soon as light travels on such null geodesies, which 
is the case in particular in general relativity and scalar- 
tensor theories of gravity. § IIII Al reviews the standard 



derivation on the propagation of a light bundle and de- 
fines the shear and convergence. In § 1111 Bl and § 1111 CI 

we apply this formalism to a Friedmann-Lemaitre space- 
time and then to a perturbed spacetime in order to get 
the final expression of the shear. 



A. General derivation 

Let us consider the evolution of a light bundle in a 
spacetime with metric g^ v following the original work of 
Ref. |3^| along the lines of Ref. 40j. The worldline of 
each geodesic can be decomposed as 



x»{\) = x»{\) + ^{\) 



(39) 



where x M (A) is a fiducial null geodesic of the bundle and 
A is an affine parameter along this geodesic. £ M is a dis- 
placement vector that labels the other geodesies of the 
bundle. The tangent vector k^ along the fiducial geodesic 
satisfies 



kfj,k^ 



0, k v WM 



0. 



(40) 



We assume that the geodesic bundle is converging at 
the observer's position O. We assume A to vanish in O 
and to increase toward the past. In O we also choose 
a quasi-orthonormal reference frame {k^, u M , n^, n% }, 
where it M is the 4-velocity of the observer and satisfies 



-1 



(41) 



and n^,ri2 are two the spacelike vectors spanning the 
plan orthogonal to the fiducial light ray, i.e. to the line 
of sight. They satisfy 



(42) 



with a,b — 1,2. From this tetrad constructed in O, we 
construct a basis at each point of the geodesic by parallel 
transporting it along the fiducial geodesic. 

As shown in Ref. |40(, £ M can be decomposed as 



(43) 



a=l,2 



and one can always choose to set £0 = 0- The propagation 
equation of £^ is obtained from the geodesic deviation 
equation as 



dA 2 ^ 



(44) 



;C = R^ a pk»k a e 

where R^ va p is the Riemann tensor and D/dX = fc^V^ 
In terms of the decomposition 143fl , this equation reads 

d 2 



(45) 



where 7Z b a = i? p va pk l 'k a n a ji n^ } is known as the optical 
tidal matrix. We have used the notation £ = £ a and 
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H£ = lZ b a £b- Decomposing the Riemann tensor in terms 
of the Ricci and the Weyl tensor, it can be rewritten as 

K» a = - l -R^k»5i + C^ a ^k a n a y b . (46) 

The linearity of the geodesic equation implies that £ 
can be related to the initial value of d£/dA through a 
linear transformation 



£(A)=2>(A)^(0). 



(47) 



Using that £(0) = (bundle converging in O) , Eq. 145|) 
gives an equation of evolution for T> h a as 



B. Background spacetime 

We now apply this general results to the case of a back- 
ground cosmological spacetime with line element 

As 2 = g^dx^dx" = a 2 \rj)g ^dx 1 * Ax" (55) 

where 77 is the conformal time and a the scale factor. We 
decompose the metric as 



g flu Ax fi Ax l/ = —A?] 2 + "fijAx l Ax^ 



(56) 



where 7^ is the metric of the constant time hypersur- 
faces. Since they are hypersurfaces of constant curvature 



dA 2 



T> = 7VD. 



(48) 



The initial conditions in for the matrix 2? are given by 



£>(0)=0 and 2?'(0)=Z, 



(49) 



with X being the 2x2 identity matrix. The matrix T> 
describes the deformation of the light bundle induce by 
the spacetime geometry. 

The direction of observation, 61, and and the angular 
position of the source Os are related to the displacement 
field by 
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s = 



g(As) 
Da(\s] 



(50) 



where A5 is the value of the affine parameter at the source 
and Da is, by definition, the angular distance of the 
source. It follows that Eq. (|47|l takes the form 



93 D^* 1 ' 



(51) 



The deformation of the shape of background galaxies is 
thus characterized by the amplification matrix 



A: 



AO, 



AO! Da (As) 
It can always be decomposed as 



A: 



1 — K — 7l 72 

72 1 - « + 7i 



(52) 



(53) 



in terms of the convergence k and of the shear 7 = 
(71772)- They can be extracted from the amplification 
matrix as 



* = i-\^A 7=2 



1 (Asa -An 
2A12 



(54) 



As was emphasized by many authors j2.1 |2J, |25j 
and demonstrated experimentally by various observa- 
tions 29], the shear can be obtained from the measure- 
ment of the shape of galaxies. 

All the definitions and derivations of this section do 
not depend on the specific form of the metric and are 
thus completely general as long as photons follow null 
geodesic and the geodesic deviation equation holds. 



7 ydzW =d X 2 + S 2 K ( X )dQ 2 



(57) 



where x is the comoving radial coordinate and dSl 2 the 
infinitesimal solid angle. The function Sk is defined by 



Sk(x) 



sin(V-K"x) sinh(> 
X, 7 



-K X ) 



K 



-K 



(58) 



respectively for K > 0, K = 0, K < 0. 

From Eq. I|57|). it is clear that Sk is the comoving an- 
gular distance so that the angular distance will be given 
by 



d a = a(i])S K (x) = a o 



Sk{x) 
1 + z 



(59) 



where ao refers to the value of the scale factor today and 
z = ao/a — 1 is the redshift. The expression of x( z ) is 
given by 



X(z) 



1 



Az 



H a J E(z) 



(60) 



where Hq is the value of the Hubble constant today and 
E(z) is defined by 



E{z) = H{z)/H = (1 + z)H/H 



(61) 



The expression for E(z) depends on the matter content 
of the universe and on the theory of gravity. It is given 
in Appendix IbI 

To solve Eq. I|48|l . we will use the fact that null 
geodesies are left unchanged by a conformal transfor- 
mation since two conformal spaces have the same causal 
structure (see e.g. Ref. |4l|). For the metric g^, Eq. I|48() 
takes the form 



dA 2 



T> = -KV. 



(62) 



The solution of this equation is trivially given by 

2> (0) (A) =S K (X)I, (63) 

where Sk is given by Eq. I|5X|I. so that 

^ (0) = I (64) 

where the subscript refers to the solutions in the back- 
ground spacetime. 
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C. Perturbed spacetime 

1. Amplification matrix 

Let us now turn to the perturbed cosmological space- 
time with a metric 



ds 2 = a 2 (77) (g^ v + h^ v ) dx^dx" . 



(65) 



Again, we work in the static conformal spacetime with 
metric g^ + h^ u . As usual, we develop the tidal matrix 
as 



T> = X» (0) + T> [1) + O (h 2 ) , 



(66) 



where involves terms of n-th order in the metric 

perturbation, and consistently we solve Eq. (|48|l order by 
order. Plugging the zeroth order solution lj6~4T) . Eq. 
reduces to 



dA 



X> (1) ^n il) (X)S K (X). 



(67) 



the solution of which is explicitly given by 



where are the Christoffel symbols. If we focus on 
scalar perturbations then the second term vanishes and 
we end up with 

T^-ab = -dab® (73) 
where the deflecting potential is defined as 



(74) 



In the case of a perturbed cosmological spacetime, it re- 
duces to 



V> - 11,0' - l:, ,0't)' 



(75) 



where 9 l is the direction of observation. This potential 
includes the effect of rotation and of gravity waves. Notes 
that in the case of pure scalar perturbation in general 
relativity and in the absence of anisotropic stress <f> = ip 
and we recover the standard result <f> = 2<p. 



X> (1) (A)=/ S K (X')S K (X- A')7?- (1) (A')dA', (68) 



so that the amplification matrix is 

A(1)(A) = r x s K (x')SK(x -y) nW{xf)dX ^ {m) 

Jo i>K(X) 

This expression is again very general and just assumes 
that photons are propagating along null geodesies. Now, 
we just need to express HS 1 ^ in terms of the perturbations 
of the metric and go back to the original (non static) 
spacetime to get the final result. 



2. Deflecting potential 

To go further, we assume that the perturbations can 
be decomposed as 



1 



h^dx^dx" = -cj)dr] 2 + Bidijdx* 



+ (-iplij + Eij) dx l dx j (70) 

where <fi an d V 7 are the two gravitational potentials in 
Newtonian gauge, Bi and describe the vector and 
tensor perturbations and satisfy 



E] = DiE\ = D,B l = 



(71) 



where Di is the covariant derivative associated to the 
spatial metric 7ij. Using that 2R pua p = h va ^0+h^ va - 
han,v(i — hvfi.ap,, we obtain that 



K§ = $h iai , afi k»Wn«<n% -~ (T a p0 g a ^Kn p b ) , 

(72) 



1 



1 d 



3. Final form 

The last step requires to go back to the original (non 
static) metric to get 



(76) 



Sk{x')Sk{x-x') 
Sk(x) 



It follows that the shear is given by by the trace of this 
matrix 

«(«.*) = l[ X SKiX l SK ( iX - X ' ) ^dx', (77) 

with A2 = d a d a . The deflecting potential is evaluated 
along the line of sight as $ = $ [Sk{x')9> x']- 

In this derivation, the source has been assumed to be 
at a given redshift, it generalizes easily to a more general 
redshift distribution, p x (x)d%, as 



K ( e ) = J Px(x)«(0)X)dX- 



(78) 



In conclusion, the shear takes the form 



k(0) 



g{ X )S K (x)A2$ [S K (X)0, X] d* (79) 



where g is defined as 



f XH , „ S K {X'-X) 
9{X) = I P(X) ^ ■ 



(80) 
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D. Comments 

Again, let us stress that this construction does not as- 
sume the validity of general relativity. It is based only 
on the validity of the geodesic deviation equation. This 
is the case in particular for all metric theory of gravity 
such as scalar-tensor theories. 

We also emphasize that using a conformal transforma- 
tion to solve the problem in a static spacetime has greatly 
simplified the derivation. 

In the late universe, when we are interested by gravita- 
tional lensing effects, the matter content is mainly domi- 
nated by pressureless matter and dark energy. It follows 
that we can define an effective density field, S c s, by the 
relation 



1 i m,0 



(a/a ) 



(81) 



where fi mj o is the value of the matter density parameter 
today. The subsecript "m" refers to the matter fields, 
including at least two components, dark matter (cdm) 
and baryonic matter (b) . If we express the 3-dimensional 
Laplacian as A2 + d 2 where d z is the derivative along the 
line of sight, then the convergence (j79(l can be expressed 
in terms of the effective density contrast as 



9(x)Sk (x) rrm d x- 



(a/ao) 



(82) 



IV. WEAK LENSING OBSERVABLES 



In the former section, we have derived the expression 
of the convergence in terms of the deflecting potential. In 
cosmology, we are interested in the statistical properties 
of the convergence. 

This section is devoted to the observable quantities to 
be compared with weak lensing data. § IIVBI describes 
the computation of the power spectrum and § IIVCI is 
de voted to the computation of various 2-point statistics, 
such as the shear and the variance of the aperture mass. 
This requires to smooth the shear field. We then describe 
our numerical implementation in § IIV Dl 



A. Power spectra 

The convergence is a function on the 2-sphere and it 
can be expanded in a 2-dimensional Fourier transform as 



(83) 



From the coefficients k(£), one can define the power spec- 
trum of the shear as 



{k(£)k*(£'))=P K (£)S^ (£-£'). 



(84) 



The shear components can be decomposed alike. It is 
convenient to use a complex notation for the shear, 7 = 
71 + *72 ■ Since the shear and the convergence derive from 
the same potential, they are not independent and it can 
be shown that the convergence can be obtained from the 
shear as 



k{£) = it (£)*/(£)/« 



(85) 



for all £ = (£1, £2) 7^ 0. The kernel function K is defined 
as 



K = Tr(£l-ej + 2ie 1 £ 2 )/£ 2 - 



(86) 



Interestingly, this implies \k\ 2 = 77* so that the shear 
and convergence have same power spectrum 

P 7 (£) = P K (£). (87) 

Note that the shear can then be obtained as 

k(0) = k + ~ J k*{0 - e')j(e')d 2 e'. 

The integration constant kq is related to any constant 
uniform mass distribution that will contribute to the con- 
vergence but not to the shear. 

Analogously, all 3-dimensional fields can be developed 
in Fourier modes as 



<t>(r,v) = / 0(fc,f?)e 



Ak-r 



dh 



(89) 



(27r) 3 / 2 

and are associated to a power spectrum 

(4>(k, v )^k', v ))^P 4 ,(k,7 1 )S^(k~k / ). (90) 

B. Shear power spectrum 



As Eqs. <|79I8U[1 show, the shear is obtained as a 
weighted projection of the deflecting potential and takes 
the general form 



«= / q(x)6 e s[S K (x)9,x]dx- 



(91) 



This implies that if ^ G ff is a homogeneous and isotropic 
Gaussian random field, so will k. Its angular correlation 
function, £ K = (k(<p + 9)n((p)), can be expressed as 



= / d x d x 'q(x)q(x') 

x(5 cS [S K (x)v,x] S eS [Sk{x){v + 0),XK92) 

The correlation function that appears in this integral can 
be computed if 5 e ff is expanded in Fourier modes l|89() to 
give 



e -ik±.ipf K (x) e -ik 3 x e -ik' ± -(<p+e)f K (x') e -ik' 3 x 



X(£eff(fc,x)<5cff(fc',x')) 



d 3 fed 3 fc' 

(2tt) 3 



(93) 
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On small angular scales (9 <§C 1), k\ ^> fc| and the power 
is mostly carried by k±. We thus approximate 

(<5 off (fc, x )<5 off (fc', x ')> - (<5cff(fc ± ,x)<5eff(fcl,X , )>'5 (1) (fc3+^). 
The integration over k% gives a factor 2tt(5 (1) (x - x') 

so 

that after integration on \' we en d up with the expression 



An interesting statistics arises when one uses a com- 
pensated filter, that is 



e«(*) = / <W(x) 



d 2 fc ± 

(2tt)2 



-P c ff(fc_L,x)e 



(94) 



We conclude, that the 2-dimensional power spectrum of 
the convergence is related to the 3-dimensional power 
spectrum of 5 e ff by 



a 1 'mat 



gix) 



a(x) 



p 



cff 



Sk{x) 



■ X 



d X - (95) 



If we filter the shear field by a window function, 
U(9' , 0) of angular radius 9 as 

«(0) = /" d 2 e'K(e')u(e',e). (96) 

then the variance of the filtered convergence is given by 
(n 2 )=2irj^ dllP K (£)(^J d'U(9',6)J (£6)de) (97) 
where Jq is a Bessel function. 



C. 2-point statistics of the shear field 

The shear field 7 has two components that can be de- 
composed in various ways. From a reference point, such 
as the center of the filter, one can define a radial and 
tangential shear as 



7r 



Re(<£*7), 7t=Im(vj*7) 



where if represents the unit vector pointing from the ref- 
erence point to the point where 7 is defined. From these 
components we can define the correlation functions 



£r(0) = (7r7r), & = (7t7t) 



(99) 



and by symmetry the correlation (7t7r) strictly vanishes. 
We can also combine this functions as 



Developing the shear in Fourier modes, one gets 



£+(0) = 



^-£p K (e)j (ee) 



<M 
2tt 



£P K (£)J 4 (£9). 



(100) 

(101) 
(102) 



Since the shear and convergence have same power spec- 
tra, we also deduce that 



df 

-7 • ->:(")= / 2^«C) 



2Ji(£9) 



CO 



(103) 



6'U{6',6)d6 = 0. 



(104) 



With such a filter, any constant mass density, related to 
the integration constant ko in Eq. (|88|l will not bias the 
statistics. It is usual to define the aperture mass as 



M, 



a p 



d 2 0'U(9',9) K (8'). 



(105) 



Interestingly, it can be expressed in terms of the tangen- 
tial shear as 



d 2 8'Q(9',9) lt (8') 



(106) 



with the filter Q defined as 

Q(9', 9) = ^ J d9' e'U{0', 9) - U(9', 9). (107) 

A widely used family of filter that satisfies the condition 
(tToHl is 



,/N 2 
1- I ~-T 



1-3 t 



(108) 



for which the aperture mass has a variance given by 



000 roo 

«,>(<?) = — / d^P K (^) 
~ Jo 



Ji{£6) 



P9 2 



(109) 



All these 2-point statistics derive from the same power 
spectrum are thus not independent. 



D. Numerical integration 



The computation of the shear power spectrum as well 
as the 2-point statistics have been implemented in a nu- 
merical code that can be used both with general relativity 
and scalar-tensor theories. It is based on the Boltzmann 
code described in Ref . , implemented to study the evo- 
lution of the background and of the linear perturbations 
for a general scalar-tensor theory specified either in Jor- 
dan or Einstein frame (see Appendix [BJ. In particular, 
this code allows the study of the magnitude-redshift re- 
lation and of the CMB angular power spectrum. It in- 
tegrates the Einstein and fluid perturbation equations 
forward in time (toward the future) and sets the initial 
conditions on the scalar field by a shooting method in or- 
der for the dark energy density parameter and the New- 
ton constant to agree with their value today (see Ref. 
for details). 

The CMB code gives access to the value of the per- 
turbation variables, among which the two gravitational 
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CMB code: 

input: f2, P(k), . . . 




Background: \{ z ) 

Perturbations: S(k,z), ip{k,z), <f>(k,z), ip(k, z), . . . 





Linear to nonlinear mapping 




FIG. 1: General description of the computation of the lensing observables. The input are generated by a CMB code that deals 
with scalar-tensor theories and extended quintessence @- Once the linear to non- linear mapping has been applied, we can 
compute the shear power spectrum and then all the lensing observables with the same normalization as used for the CMB 
observables. 



potential 4> and ip and the density perturbation, as a func- 
tion of the wave numbers k and the redshift z. There- 
fore we do not need, as is usually done in the Newtonian 
regime, to decompose the density perturbation in an ini- 
tial random field and a growth factor. In such a case, 
the lensing observables depend on the normalization and 
shape of the transfer function that is on as and T. 
In our approach, all observables are CMB normalized at 
£=111, according to Ref. > an d we need not introduce 
these parameters. The CMB code also gives the value of 
from which we can derive the angular distances. 
This implies that the lensing plug-in code (see Fig. 
we developed, does not need to solve any evolution equa- 
tions. It only integrates the deflecting potential along 
the lines of sight accounting for the source distribution, 
according to the two-dimensional projection described in 
Sec. IIV Al and deal with the linear to non-linear mapping 



rsees iym^ . 

We are interested in angular scales ranging from 1 ar- 
cmin to 2 degrees and we considers multipoles in the 
range [£ m in = 2,^ max = 22000]. We thus have to consider 
comoving wave numbers in the range [fc m i n , k max ] with 

fcmin ~ 4nn/x(^max) ~ W^h Mpc" 1 (110) 

and we set the cut-off A: max ~ £max/x( z miu) to 10 h 
Mpc . Since some of the scales are in the non-linear 
regime, we carry out the linear to non-linear mapping 
applying the procedure of Ref. [44| . It has to be stressed 
that these mappings have been calibrated on numerical 
simulations assuming that gravity was described by gen- 
eral relativity and ACDM cosmological models. The va- 
lidity of this hypothesis will call for further checking, but 
it is to be expected that the scalar field does not affect too 
much the clustering of the matter during the non-linear 
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phase and that its main effect is through its contribution 
to the expansion of the universe, as we will discuss later. 

We assume a redshift sources distribution parameter- 
ized as 

with (z s ,a,f3) — (0.8,2,1.5) consistently with a limited 
magnitude Iab = 24.5 (for details concerning this choice, 
see Ref. 0). 

E. Non-linear power spectrum 

The computation of weak lensing observables requires 
to determine the power spectrum in the non-linear 
regime. 

For ACDM models, various mappings have been pro- 
posed in the literature. Assuming stable clustering, it 
was argued [|(| 0] that the effects non-linear evolution 
can be described by a mapping between the linear and 
non-linear power spectra involving a universal function, 
f n l. Introducing A as 

A 2 (fc) = 4irk 3 P(k), (111) 

the non-linear power spectrum is obtained by 

A 1 2 ll (fc nl ) = /„i[A 2 (fc)] (112) 

where the wavenumber, fc n i, is related to the linear 
wavenumber, k, by 

k 3 = [i+Ai l (k nl )r i k 3 nl . (us) 

The function f n \ is determined by A^-body simulations 
and it depends on the value of the cosmological parame- 
ters. It has also been shown 0] that at large values of 
its argument this function behaves as 

/nl(x) ~ 3;r3/2 ' (114) 

D being the linear growth factor [see Eq. (|118fl ] . Because 
of this simpler asymptotic analytic form, we will use the 
relations (|111I113J) for general analytic arguments. Nu- 
merically, we have however implemented the more real- 
istic mapping described in Ref. 44] . Fig. [21 depicts the 
linear and nonlinear power spectra of the density pertur- 
bation computed for the fiducial ACDM model. 

As we have already stressed, all the mappings have to 
be calibrated on A'-body simulations. No full ./V-body 
simulations for quintessence models and have been per- 
formed so far (see however Ref. 0] where A^-Body sim- 
ulations with a modified expansion rate to take into ac- 
count quintessence have been investigated). However, 
as it has been shown in Ref. |4!j, the shape of the lin- 
ear power spectrum for sub-Hubble modes, but not its 



absolute amplitude, of quintessence models is very sim- 
ilar to the one of a pure ACDM. Hence we assume that 
these mapping also apply to quintessence models. In par- 
ticular, it was argued [5(| that the mapping I|11H113[) 
was reasonably accurate for effective quintessence mod- 
els, at least at low redshift. In particular, we assume that 
the non-linear regime always reaches a stable clustering 
regime. Notice that various mappings lead to results that 
agrees only at a 5-10% level (see e.g. Ref. 0). 

As far as scalar-tensor theories are concerned, as for 
quintessence, there is no A^-body simulations to calibrate 
the mapping. We have to assume that the mapping 0] 
calibrated on pure ACDM still hold. It can be argued 
that we do not expect the change of the theory of gravity 
to drastically affect this mapping as long as strong field 
effects, such as spontaneous scalarization |20j . appear. 
Even though, we emphasize that the mapping procedure 
has to be adapted in the case of scalar-tensor theories, as 
it will be discussed in see § IVI B 21 



100000 




0.1 1 — — — 1 

0.01 0.1 1 10 
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FIG. 2: Matter power spectrum P m (k) for a flat fiducial 
ACDM model defined by Q A = 0.7, Q b h 2 = 0.019, h = 0.72. 
Linear (solid) and non-linear (dashed) regimes are presented. 



V. WEAK LENSING IN GENERAL 
RELATIVITY 

A. Generalities 

When gravity is described by general relativity, then at 
late time matter dominates and the anisotropic stresses 
are negligible so that the two gravitational potentials are 
equal 

</> = V (115) 

consequently the deflecting potential is simply given by 

$ = 20. (116) 

On sub-Hubble scales, the gravitational potential is re- 
lated to the matter density perturbation by the Poisson 
equation which implies that 

Sett = Sm: (H7) 



12 



10"' 
10" E 
1CT 
10" 1C 

10-" 
io- 1£ 

10" 13 
10" 14 



ACDM (lin) 

ACDM (nl) 




ACDM (lin) 

ACDM (nl) 



10 

) (arcmin) 



FIG. 3: (left) Convergence power spectrum P K (l) as a function of the multipole £ and the 2-point statistics of the shear field: 
(middle) the aperture mass variance and (right) the shear variance in the case of the fiducial ACDM model 1123H . 



where S m — Sp m /p m is the matter density contrast. Let 
us stress that on large scales, one must take into account 
the radiation anisotropic stress which will induce a de- 
parture from this equality (see Fig. I13[) . Starting from 
an initial time, a^, where the modes of interest are sub- 
Hubble, one can decompose the density field as 



S m (k,a) = L»(a)(5(fc,a l ) = D(a)Si 



(118) 



where D(a) is the growth factor. It follows that 
P e g(k,a) = D 2 (a)Ps{k). Its equation of evolution, in 
the Newtonian regime, is given by 

D + HD - knGp m a 2 D = 0. (119) 

It can be rewritten by using the redshift as variable as 



D" 



- 7-^) D ' - |fim,o(l + z)D = 0, (120) 



setting ao = 1 and where a prime refers to a derivative 
with respect to z. This equation has two solutions, a 
decaying mode, D oc H, and a growing mode 



A tin 



da 



-\i3 " 



\aE(a)\ 



(121) 



The convergence power spectrum takes the simplified 
form 



Sk{x) 



(122) 



Figure |2| depicts the convergence power spectrum, the 
aperture mass and the shear variance for a ACDM model. 
This model will be our reference model and we will try 
to quantify the deviation from its prediction on various 
models. 



C. Quintessence: effect of the potential 

As a first generalization to the previous ACDM model 
we consider a class of quintessence models (QCDM) with 
runaway potentials of the form 



V(ip) = M 4 (p- 



(124) 



In the slow-roll regime, the quintessence field acts as a re- 
pulsive matter component that replaces the cosmological 
constant. Gravity is still described by general relativity 
and the main effects on the lensing observables arise from 
the modification of the Friedman equations, and thus of 
the angular distances and of the growth factor of the 
density field. 

To estimate the amplitude of the effects, let us follow 
Ref. [z| and assume that the sources are located at a 
redshift z s so that the source distribution is simply given 

byp(x) = 6(x~Xa) and thusg(x) = S K (Xs-x)/ 'Sk(Xs)- 
The function 



W(x,Xs) 



Sk(x)Sk(Xs - x) 



(125) 



Sk(Xs)ci 

is peaked around Xs/2 so that it can be approximated by 



B. ACDM 

The power spectrum l|122|l depends on the cosmologi- 
cal parameters through the growth function and the an- 
gular distances and on the normalization of the power 
spectrum. As emphasized above, we do not perform such 
a splitting and use the complete expression 1)95(1 so that 
our results depend on the cosmological parameters and 
the primordial spectrum are CMB-normalized. 

As a reference model, we choose a flat ACDM with 

O\,o = 0.7, il m .o = 0.3, a,. = 0.037, h = 0.72. (123) 



W 2 ~ W 2 /2 S(x - X1/2) 



with X1/2 - Xs/ 2 and 



Wt/2 = 



W 2 (x,Xa)d X . 



(126) 



(127) 



Plugging this approximation into Eq. (|95|l allows to esti- 
mate the shear power spectrum as 



Wl/2 



Sk(xi/2 



P 



Sk(Xi/2 



IX1/2 



(128) 
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FIG. 4: Comparison between the fiducial ACDM model l|123|l and various quintessence models, (left) the growth factor 
normalized to its value in an Einstein-de Sitter universe, D(a)/a, as a function of the scale factor a, normalized at high redshift; 
(middle) Relative deviation from ACDM on the comoving radial distance, x( z )i antI (right) on the geometrical factor, q(z). In 
all plots, we consider three QCDM models with potential 1124H with m = 6,8, 11 from top to bottom. 



Assuming that the initial power spectrum takes the form 



Ps(k) = Ak n , 



(129) 



we get that in the linear regime P in = D 2 Ak n and P n i 
a 3^3(3+n)/2 SQ that Eqg (lining i mp iy that 



k 



(ejAa*) l/{5+n) 



(130) 



As long as we are considering modes in the linear regime, 
the effect of the quintessence field is mainly encoded in 
the growth function. We get 



Pi 



W. 



Q " 

1/2 



lin 



s Q Ax Q 1/2 ), 

'D? /2 /D%\ 




2+n 



pA 



(131) 



where D\i<i and Dq are the values of the growth fac- 
tor at % = X1/2 and Xo respectively and Pq i the 
value of the matter power spectrum today, evaluated at 
k = 1/iMpc -1 . 

In the non-linear regime, the matter power spectra are 
not distorted in the same way because the same value of 
fc n i does not correspond to the same k. It follows that 



pQ 

K 
P A 



111 




3/(5+n) 



(132) 



if the spectral index is the same in both models. The rela- 
tions (|131fl and 1|132|) show that on small scales the shape 
of the matter power spectrum is modified in quintessence 
models. 

As can be seen on Fig. matter perturbations grow 
more slowly than in a ACDM model. The CMB nor- 
malization implies that at lower redshift the amplitude 
of density perturbations are smaller in quintessence mod- 
els. Thus, a given scale enters the non-linear regime later 



TABLE II: Quintessence models (QCDM): background ef- 
fects. (Upper table) Maximum relative deviation on the co- 
moving radial distance, x( 2 )> geometrical factor, q(z), and 
linear growth factor, D(z), from the fiducial ACDM model 
for three models with inverse power law potential, Eq. 111241 . 
(Lower table) Quantities used to estimate the amplitude of the 
effects on the weak lensing observable, according to Eq. 11311 . 
The amplitude of the matter power spectrum Po is evaluated 
at k = 1 /iMpc -1 . For the ACDM model zy 2 ^ 0.51, for all 
the quintessence models 2]/ 2 — 0.48. 



m 


Ax(z) (%) 


Aq(z) (%) AD{ 


z) (%r 


6 


-17 at (z = 3.4) 


-20 at (z = 1.4) 


-42 


8 


-19 at (z = 3.9) 


-23 at (z = 1.6) 


-51 


11 


-21 at (z = 4.7) 


-26 at (z = 1.7) 


-61 


m 


Xi/a/Xi/2 


/ 2 /(Wi 2 / 2 ) A D 1/2 /D 


Po/P A 


6 


0.853 


0.611 1.052 


0.275 


8 


0.844 


0.580 1.056 


0.189 


11 


0.835 


0.549 1.060 


0.113 



"Evaluated at z = 0. 
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FIG. 5: The ratio of the three dimensional power spectrum 
of the matter perturbation P m (k) of two quintessence models 
to the reference ACDM model 11230 for m = 6 and m = 11. 



than in a ACDM. This time delay accounts for the sharp 
modification of the spectrum at small wavenumbers. Ta- 
blelTTI gives the order of magnitude of the expected effects 
on the background quantities, the growth factor and the 
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TABLE III: Quintessence models (QCDM): absolute values of the convergence power spectrum, P K (£), aperture mass variance, 
(Ma P (@))i an< i shear variance, (7 2 (#)) at two angular scales respectively in the linear and non-linear regimes. For the small scales 
we give both the value in the linear and non-linear regime (values within parenthesis). Each model is labelled by QCDMm, 
where m defines the inverse power law potential, Eq. 112411 . 
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"The multipoles I = 180, 7200 correspond to the angle = 2°, 3'. 
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FIG. 6: (left) Convergence power spectrum P K (£) as a function of the multipole £ and the 2-point statistics of the shear field: 
(middle) the aperture mass variance and (right) the shear variance for the reference ACDM model l|123|l (solid/long-dashed 
lines) and two quintessence models (m = 6, short-dashed/dotted line; m = 11, long-dashed-dotted/short-dashed-dotted line). 



lcnsing observables. For instance let us concentrate of 
a quintessence model with m — 6. Table ITTT1 shows that 
the convergence power spectrum is 80% to 85% smaller in 
quintessence than in ACDM. Now, assume we normalize 
the power spectra on linear scales at z = 0, it will have 
implied to multiply the linear power spectrum by ~ 7 
and the non-linear one by - V7. Thus, we would have 
found that the non-linear regime will differ roughly by 
50% from the ACDM. This conclusions are similar to the 
ones of Ref. and show that the change of spectrum 
between the linear and non-linear part of the spectrum 
sets strong constraints on the time evolution of the dark 
energy. 



VI. WEAK LENSING IN SCALAR-TENSOR 
THEORIES 

We now turn to the case where gravity is not described 
by general relativity but by a scalar-tensor theory. 



A. Newtonian regime 

Before we discuss explicit models, we can try to eval- 
uate and discuss the expected effects on lensing observ- 
ables. For that purpose, let us first look at the pertur- 
bation equations in the Newtonian regime. We consider 
modes with wavelengths smaller than the Hubble length 
and also assume that the scalar field is light. 



In the matter era and on sub-Hubble scales, the equa- 
tion (|B6p of Appendix 151 reduces to 

^ - <P = (133) 

so that non-minimal coupling induces an extra- 
contribution to the anisotropic stress, while Eq. (|B8|) re- 
duces to a generalized Poisson equation 

FA(f) = A-KG*pa 2 8 m - -^-ASip. (134) 

The evolution equation l|B4(l of the matter field takes the 
form 

5 m = -AV + 3ip (135) 

and the Euler equation (|B5|I to 

V + HV = -(j). (136) 

This set has to be completed by the Klein- Gordon equa- 
tion for the scalar field evolution 

(A - U w a 2 ) dip = F V A(<P - 2yj). (137) 

B. Effects and amplitudes 

Various effects are expected. As in the case of 
quintessence, the Friedmann equation is modified so that 
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the geometrical function q as well as the growth of density 
fluctuations are modified. 

Besides, there are two specific effects which do not ap- 
pear in pure quintessence models due to the modification 
of Einstein equations: 

1. The Poisson equation does not take its standard 
Newtonian form so that 



(138) 



2. The two gravitational potentials are not equal any- 
more so that 



$ ^ 20. 



(139) 



By extension to the post-Newtonian formalism, we 
can define the parameter ^y PPN in the cosmological 
context by 



V = 7 PPN 
so that $= [l + j FPN (z)](j>. 

1. Linear regime 



(140) 



The density evolution follows an equation that is sim- 
ilar to the pure Newtonian case but where the gravita- 
tional constant has to be replaced by its time dependent 
value 



HS m - 47rG C avpa $m = 0. 



(146) 



If we decompose the density field as in Eq. (|118fl . we 
obtain that the effective growth factor is 



D eS ^^[l + r] ] 2 D(a) 



(147) 



where D is solution of Eq. I|14tj|) ■ Written in terms of the 
redshift, it takes the form 



D" 



IT 



1 



d' --<),„ od f n i). 

(148) 

where fi m o is defined by Eq. l|A"l9)) . It follows that the 
shear power spectrum takes the form 



9H 4 n 



n 



m,0 



g(x) 

a(x) 



xD 2 ( x )Ps 



[i + vY 



Sk(xY 



X 



d X . (149) 



The two equations of evolution (|135I136[1 can be com- 
bined to get the evolution of the density field 

S m + 2HS m -A<t> = 3a 2 (jm- 2 ^) ' . (141) 



It follows from this expression that we can estimate P K 
by solving the background equations solely. In particular, 
because of the form i|144|) of the Poisson equatiom we do 
not expect effects as the ones described in Ref. [4(| on 
small scales (see however Fig. El- 



When Utptp is much smaller than the wavelength of the 
modes we consider then we can combine Eq. (I133|l and 
Eq. I|137l) to get the fluctuation of the scalar field as 



dip 



FF,„ 



F + 2F$ 



In particular, this implies that 

F 2 

7 PPN (z) - 1 = v 



F + 2F 2 



(142) 



(143) 



2. Remark on the linear to non-linear mapping 

The linear to non-linear mapping procedure has to be 
extended when working in scalar-tensor theory, in par- 
ticular because we will have to deal with the scalar field 
perturbations. As Eq. H95|l shows, we must determine 
P G "g. For that purpose, we use the definition llHT|) to de- 
fine S c f[ and we decompose it as 



(150) 



similarly, as obtained in § IIIBI in a non-cosmological con- 
text. Inserting Eq. 1142|) in the generalized Poisson equa- 
tion l|134fl . one gets 



A4> ~ AirG c&v pa 2 8 n 



(144) 



where G cav is defined by Eq. l(2!)|) and is indeed time 
dependent. Using the Eq. I]8ip. this equation implies 
that 



Fo 
F 



(145) 



where rj(<p) = 2F 2 /(2F + 3F 2 ). Today r)Qpo) is expected 
to be smaller than a few 10 -3 . 



where 5x contains the contribution of the scalar field 
perturbation and its derivative. 

Assuming that the scalar field does not cluster, as it is 
the case in quintessence, Sx will not enter the non-linear 
regime so that we can assume that 



5 X = 6 1 ™ 



(151) 



In the Newtonian linear regime, Eq. I|81|) implies that 



rim _ 
°cff — 



a/a 



(152) 



The output of our CMB code gives access to $ lm and <5jj" 
from which we can deduce <5^ and thus S 1 ^ 1 . With this 
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ansatz, we get the effective perturbation in the non-linear 
regime 



jrnl c-nl , jrlin 
<W = V + °X • 



(153) 



It follows that 



P^(k, z) = P^(k, z)+P^(k, z) + 2y/p$(k, z)P^(k, z). 

(154) 

When we are deeply in the Newtonian regime, but still 
in the linear regime, Eq. (|H5j> implies that S x -> {F /F- 
l)S m so that 



P^(k,z) 



neglecting the contribution in 77. In fact, notice that the 
ratio Pig I P 1 ™ evaluated today is always grater than one, 
because of the non vanishing anisotropic stress of the 
radiation and non- minimally coupled scalar fields. In the 
non-linear regime, 



Fo 
F 



- 1 



so that at z = 0, P^(k) — > P'^(k), as expected from our 
assumption that the scalar field does not cluster. 

Hence, our ansatz amounts for an interpolation be- 
tween the super-Hubble regime where the contribution 
of the scalar field perturbation cannot be neglected and 
the Newtonian regime where the effective Poisson equa- 
tion (|144fl hold. Note that it assumes that the scalar 
field does not enter the non-linear regime. In fact the 
mapping procedure of Ref. uses a halo model and 
it was shown |5l| that in a spherical collapse, the scalar 
field does indeed not cluster (see also Ref. [E3,IE3])- 
we can hope this interpolation to be justified but, at the 
time being, we have no possibility to check it. 



3. Effect on the shear power spectrum 

Following the approximation of Section IV CI we can 
compare two models which differ only by the theory of 
gravity (i.e. same quintessence potential or ACDM mod- 
els). 

In the linear regime, we get 

iV ^xf/p y^/ fo 

' {S S K T (X?J 2 )) \F(z 1/2 ) 



P 



ST 



pRG 



Iin 



'1/2' 

while in the non-linear regime we obtain 

2 / ~ r>n t \ 2 



pRG 



(155) 
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(.156) 



The main contribution to the modification of the shear 
power spectrum, and consequently on the 2-points statis- 
tics, is expected to arise from the evolution of the matter 
power spectrum. Because of the normalization to the 
CMB, i.e. at high redshift, its amplitude evaluated to- 
day accounts for integrated effect over a wide redshift 
range. Therefore it leaves trace on the shear power spec- 
trum, even if the computation of the shear power spec- 
trum involves an integration just up to the source redshift 
z s . Instead, as for quantities depending on the redshift 
range [0, z s ], eventually evaluated at 21/2, their contri- 
bution can become relevant only if Z1/2 1S sufficiently 
high. In particular, for Z1/2 — 0.5 the deviation from 
general relativity described by F /F 1 / 2 is negligible (see 
Fig. [7| hence the peculiar effects of scalar-tensor theories 
are ultimately encoded in the amplitude of the matter 
power spectrum evaluated today. However notice that, 
accounting for sources at higher redshift, where the de- 
viations from general relativity could become significant, 
and looking at lenses at low redshift, one should detect 
deviation of order 10% on weak-lensing observables, pe- 
culiar of a scalar-tensor theory. 



VII. INVESTIGATION OF TWO EXPLICIT 
MODELS 

We consider two test models to discuss in more de- 
tails the amplitude of scalar-tensor theories on lensing 
observations. Indeed, scalar-tensor theories have two free 
functions that need to be specified, which open a pa- 
rameter space much larger than a standard ACDM. Any 
quintessence model is specified by its potential U(ip). To 
embed it in scalar-tensor theories, one has two possibili- 
ties: (1) assume that the potential in Jordan frame, U(ip) 
is the quintessence potential and choose a coupling func- 
tion F{<p) or (2) assume that the potential in Einstein 
frame, V(ip*), is the quintessence potential and choose a 
coupling function a(y>*). The second possibility is prob- 
ably more secure from a theoretical point of view since it 
amounts to specify the property of the true spin-0 degree 
of freedom of the theory. 

We will however investigate the two possibilities. We 
first consider in S I V 11 Al the case of a non- minimally cou- 
pled scalar field and in § IVII Bl the case of scalar-tensor 
models that incorporate attraction toward general rela- 
tivity. 



A. Non-minimally coupled case 

We consider the simplest model we can think of for 
which the coupling of the scalar field to gravity is de- 
scribed by the function 

F( V ) = l + & 2 . 



The constraint l|25|) implies that 

|f Pol < 0.55 x 10" 3 . 



(157) 
(158) 
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FIG. 7: Deviation from general relativity parameterized by 
7ppn(z) — 1 for several non-minimally coupled models with 
Ratra- Peebles potential, Eq. 1)1240 . as a function of red- 
shift z. Solid line (long dashed line): non-minimal coupling, 
£ = +0.001 and m = 6 (m = 11) [potential defined in Jor- 
dan frame]. Short dashed line (dotted line): exponential 
Damour-Nordvedt type coupling, with = 4, B = 0.5 and 
m — 6 (m = 11) [potential defined in Einstein frame]. Hori- 
zontal lines: level of the upper bounds measured by gravita- 
tional experiment on Solar System scales [z ~ 0). 



Generically 0, tpo 
implies that 



0(1) today so that this constraints 



|f| < 0.5 x 10 



-3 



(159) 



In this class of models, all deviations from general rela- 
tivity will scale like £ so that we expect 5P K /P K ~ 0(£), 
that is very small effects. 



1. £ACDM: effect of the coupling 

To quantify the effect of this coupling, we first assume 
that tp has no potential and that there exists a cosmolog- 
ical constant. 

At this point, we should stress that there is no unique 
way of defining a cosmological constant in scalar-tensor 
theories. Either one adds this constant in Jordan frame, 
in which case the associated energy density is constant. 
But, from Eq. @, this induces a potential for the true 
spin-0 degrees of freedom of the theory. On the other 
hand, one could argue that imposing a constant poten- 
tial in Einstein frame lets the scalar field massless hence 
generalizing a constant potential. But, it will not corre- 
spond to a constant energy density in the physical (Jor- 
dan) frame. Let us also mention that adding a cosmolog- 
ical constant in Jordan frame modifies the potential as 
U — ► U+ A, corresponding to the change V — > V + AA i /2. 
It follows that the spin-0 degree of freedom can remain 
massless only if 

V"fo>.) + 2AA 4 (<^) [4a 2 + 0\ = 0. 

We have considered various models with U{<p) = A 
and with the constraint Q\ = 0.7 for various value of 
£, typically ranging from —0.005 to 0.1. None of them 



exhibit any departure from the standard ACDM reference 
model <|123f) . as expected from our general arguments. 



2. Extended quintessence 

We now take into account both effects of a non- 
minimal coupling and of a runaway potential in Jor- 
dan frame. We consider models with m = 6, 11 and 
£ = 10" 2 , 10" 3 , 10" 4 , -5 x 10~ 3 . We will refer to these 
models by the label EQm — £n, where m stands for the 
exponent of the runaway potential and n for the value of 
the coupling (n = 1, ... ,4 for £ = 10~ 2 ... - 5 x 10~ 3 , 
respectively) . 

As expected, the main contribution to the deviation 
on the shear power spectrum and consequently on the 2- 
points statistics come from the amplitude of the matter 
power spectrum evaluated today. The relative differences 
occurring in the (comoving) angular distance, in the win- 
dow function and in the growth factor at z — z±/2 — 0.5 
are negligible. Notably, the evolution of the coupling 
function F(tp) is negligible. Anyway, the scaling in £ is 
obvious (see Tables EIIIU and EIIH3- 

As can be shown on Fig.|S|we obtain an effect of order 
20% for £ = 0.01, a value excluded by Solar System con- 
straints, see Fig. |7| For £ = 10~ 3 , our results agree with 
those of Ref. [3(j- They argue that, taking into account 
the different normalization of the scalar field compared 
to our work, 

SP K /P K - -w&l 

hence leading to a ~ 2% effect on the shear power spec- 
trum once the constraints (|158I159|I are taken into ac- 
count. In general, since the effect scales like £, we can 
conclude that for values compatible with Solar System 
constraints, the amplitude of the effects will be smaller 
than a percent. Such a low value makes these models, 
in practice, indistinguishable from their corresponding 
quintessence models. 



B. Attraction toward GR 

1. Mechanism 

Another interesting class of models is the one in which 
the scalar-tensor theory is attracted toward general rela- 
tivity today. This feature is better described in Einstein 
frame. The Klein Gordon equation (cf. Eq. IA9|1 can be 
rewritten in terms of the new variable 



P 



In a* 



(160) 



as 



7^— ip* + [l-w + 2v{y>*)]<p lf 



-a(^*)(l - 3w) - v{(p*)^^- (161) 
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FIG. 8: Comparison between an extended quintessence model with £ = 10 -2 and m = 6 and a quintessence model with same 
potential, (left) convergence power spectrum ([P k ]eq-£/[-Pk]qcdm — 1), (middle) aperture mass variance and (right) and shear 
variance. Solid (long dashed): linear (non-linear) regime for Ratra-Peebles potential with m = 6. Short dashed (dotted): linear 
(non-linear) regime for Ratra-Peebles potential with m = 11. 



where a dash temporarily refers to a derivative with re- 
spect to p, w refers to the equation of state of all matter 
fields but <y9* and where we have introduced the reduced 
potential 



The coupling function A(tp*) is decomposed as 



A(<p, 



(162) 



(163) 



so that a = da/dtp*. 

The mechanism of attraction is then well illustrated by 
the original model [l|| in which the scalar field potential 
is fiat and where a(cp*) is quadratic. Setting a(ip*) — 
/3<f%/2, the Klein-Gordon equation takes the form of the 
equation of motion of a particle with velocity-dependent 
inertial mass, m(ip^) — 2/(3 — V* ) j subject to a damping 
force —(1 — w)ip'^ in a potential, (1 — 3w)a(tp*), 



- ^': + (l-w)^=-f3(l-3w)^. (164) 

The positivity of the energy implies that m(<£>*) > 0. 
During radiation era w = 1/3 and the field is decoupled 
from the potential and will tend to a constant value, 
<p# t R, whatever its initial velocity. During the matter 
era, the evolution of tp* is the one of a damped oscilla- 
tor that starts with a vanishing initial velocity, tp* will 
thus moved toward the minimum of a(<p*) where a = 0, 
if /3 > 0. That is the scalar-tensor theory will become 
infinitely close to general relativity. 



2. Massless scalar field models 

We first investigate a class of models with a vanishing 
potential and a quadratic coupling defined in the Einstein 
frame, as in the original Ref. 15] . The coupling function 
takes the form 



a(tp*) = a m + i/3(</?« - ip m ) 2 



(165) 



where tp m is the value of the field at minimum. The 
functions a and f3 defined in Eqns. I|13|l and H14[) are 
thus given by 

a(tp*) = /3((p* - tp m ), (3(tp*) = f3. (166) 

The coupling i|165[) can be rewritten as 

a(ip*) =a + a (ip« - ipQ + -0(<p* - (p ) 2 (167) 

where (po is the value of tp* today. Without loss of gener- 
ality, a redefinition of units allows to reduce to this form 



= (168) 
The post-Newtonian constraints (|23I24(1 imply that 



a 2 < KT 



-8.5 x 1CT 4 < al(l + Po) < 1.5 x 10" 



(169) 

A detailed analysis [16j of the primordial nucleosynthesis 
in this particular model also give constraints on (a 0j /3)- 
During the matter era, assuming tp'J 1 -C 3, the Klein- 
Gordon equation simplifies to 



'-y'l + & + pep* = 0. 



(170) 



According to the value of /?, it has two different kinds of 
solutions. For (3 < 3/8, that is for a small curvature of 
the coupling, 

p.(z„) = A+(1 + z^ 1+r ^ 4 + + z*) 3 ^/ 4 (171) 



where we have set r = yl — 8/3/3. The solution com- 
patible with the bound l|169[) is thus given by 



.(**) = ^(1 + z*) 3(1_r)/4 . 



(172) 



For j3 > 3/8, we have damped oscillations around the 
minimum of a(<p*) 



¥>•(«.) = (1 + z*) 3/4 {a< 
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In this situation, G cav may be oscillating during the mat- 
ter era. In the radiation era, the solution for the scalar 
field can be shown 15! to behave as 



<P 



* = ±V3arctanh\/l — A(l + z)~ 



B. 



(174) 
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FIG. 9: Evolution of the scalar field in Einstein-frame in a 
model of massless dilaton with a quadratic coupling (qo = 
10" 4 , P = 10 2 ). We compare the case of a vanishing potential 
(solid) and a constant potential (dash). The solutions differ 
only recently. The period of the oscillation is compatible with 
the analytic result, as can be seen on the left plot where <p* 
has been divided by the damping dactor (1 + z*) 3 / 4 . 

In figure we show that the numerically computed 
evolution of the scalar field fits very well the previous 
analytic solutions. Interestingly, the models considered 
in Fig. |21 are compatible with nucleosynthesis and Solar 
System bounds. We can check that the oscillation of the 
scalar field do not leave any imprint on both matter and 
CMB angular power spectra. 
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FIG. 10: Relative deviation (P k /[P k ]scdm - 1) from a SCDM 
model (A = 0) on the convergence power spectrum of the 
massless scalar fields models, with quadratic coupling with 
(ao,f3) = (10 _4 ,0.1), as a function of the multipole I. Linear 
(solid) and non-linear (dashed) regime are presented. The 
field is attracted rapidly toward the minimum of the coupling 
function so that the model does not differ from a SCDM by 
more than 0.1%. 

Figure lTUI depicts the shear power spectrum for a model 
with a massless dilaton with quadratic coupling 1)165(1 
with (aQ,/3) = (10 _4 ,0.1). The scalar field is attracted 
toward the minimum of the coupling very efficiently so 
that the model can hardly be differentiated from the anal- 
ogous model in general relativity. We emphasize this 
class of models, which have also been constrained by 
BBN 16], do not account for the late acceleration of the 



universe. Thus, we have compared them to a standard 
cold dark matter model with A = 0. 

Note that adding a constant potential, that is a cosmo- 
logical constant defined in Einstein frame, lets the dilaton 
massless and does not affect its dynamics apart from at 
very small redshift (see Fig.[5J). 

As for non-minimal quadratic model, this class of 
model will not leave any significant signature on lensing 
observables. 



3. Extended quintessence 

In the previous models, the scalar field accounts for a 
new interaction but not for the acceleration of the uni- 
verse that was driven by a pure cosmological constant. 

In the context of quintessence, extended quintessence 
models have been widely considered. In particular, it was 
realized that scalar-tensor quintessence models with 
attraction toward general relativity can be constructed. 
This is the case for instance when quintessence is con- 
structed as a runaway dilaton |l8l Il9|. In this case, the 
coupling function typically reduces to 



-Be-to' 



and the potential takes the form 



(175) 



(176) 



During radiation era, the coupling is not effective and 
the scalar field evolution will get its standard attractor 
solution. In the matter era, the field will start by slow 
rolling. If we assume it is always slow-rolling and that it 
explains the acceleration of the universe today, then the 
Klein-Gordon equation reduces to 



Be 



whose solution is well approximated by 



It follows that 



Q'O 



1 - 0a o ln(l + z») 



(177) 



(178) 



(179) 



where a — —Be~too anc [ p Q — —f3a . Besides, as in 
quintessence, the value of ipo, is obtained from the con- 
straint on f^A- Note also that no clear bounds from BBN 
but the general ones l|31l) have been inferred. Another 
constraint on (p' arises from the bounds on the time vari- 
ation of the gravitational constant. 

To illustrate this scenario, we consider four flat 
quintessence models with fi^o = 0.7, runaway potential 
(|T7oT> with m = 6. 11 and 

1. B = 0.5, (3 = 4, labelled by DN4; 

2. B = 0.1, p = 4, labelled by DN5. 
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FIG. 11: Ratio D(a)/a as a function of the scale factor a, normalized at high redshift, for extended quintessence models with 
quadratic (£ = +0.01, left panel) and exponential coupling (/3 = 4, B = 0.1, right panel) and for the corresponding minimally 
coupled models. Notice that the potential are defined in different frame. 
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FIG. 12: Comparison between an runaway dilaton model with /3 = 4, B = 0.1 and m = 6, 11 and a quintessence model with 
same potential, (left) convergence power spectrum {[Pk]eq~^/[Pk]qcdm — 1), (middle) aperture mass variance and (right) and 
shear variance. Solid (long dashed): linear (non-linear) regime for Ratra-Peebles potential with m = 6. Short dashed (dotted): 
linear (non- linear) regime for Ratra-Peebles potential with m = 11. 



The two choices of the coupling parameters are chosen 
to satisfy the "old" Solar System and "new" Cassini 
constraints, respectively. These models are labelled by 
EQm — DNn, where m defines the inverse power law po- 
tential and n = 4, 5. We have to stress that the deviations 
on BBN constraints for such a class of models need to be 
studied. 

We compare each model with its related quintessence 
model. Figure depicts the effect of the coupling on 
the growth of density perturbations. As can be shown 
on Fig. 1121 this implies effects of more than 10% on the 
shear power spectrum and thus on all the 2-point statis- 
tics. On the scales of interest, we can safely use the 
relations quoted in S IV1B3I Tables IVlTBll and IV11B31 
show that the amplitude of the matter power spectrum 
evaluated today is the main source of deviation on P K . 
Indeed this amplitude takes into account the whole his- 
tory of the modes since the CMB, and in particular in 
eras where large deviations from general relativity are 
now possible, see Fig. Interestingly, in the shallow 
universe, deviations from general relativity are small and 
effect on the background quantity around Z]y 2 — 0.5 are 
of the order of a couple of percents. The comparison of 
the number presented in tables IVII B 31 and IVII B 3l with 
the estimate (|155f) in the linear regime agrees within a 
factor of a few percent. 

For instance, the model labelled by EQ11 — DN5 fits 



the Solar System (Cassini) constraints and lead to a de- 
viation of order 10% from the corresponding minimally 
coupled (QCDM) model on the shear power spectrum. 

Besides, there are other signatures of these models de- 
scribed in Fig. El First, the Poisson equation differs from 
its standard expression since the scalar field, which is a 
dark component contributes to the total energy density. 
This implies that 5 c s 7^ S m . Indeed on small scales we 
have found that they are proportional [see Eq. (|144fl ] and 
that the proportionality factor is ^-independent. This is 
not the case on large scales anymore. This implies that 
the shape of P m and P e ff will differ on small k. Lensing 
is sensitive to P ff while galaxy catalog may determine 
P m . This effect is of order 10% and can be hoped to be 
detected, as first pointed out in Rcf. £4| . Second, in this 
model the CMB angular power spectrum will be modi- 
fied: the integrated Sachs- Wolfe effect will be amplified 
and the amplitude of the secondary peaks will be smaller. 

Note however than on large angular scales, P e ff ^ Pm 
even for a pure ACDM model (Fig. I13|) . The reason is 
that on these scales, one cannot neglect the anisotropic 
stress of the radiation (see Eg. IB6(I . Even if the contri- 
bution of the scalar field involves a deviation from the 
usual Poisson equation, it is not clear that this effect is 
not blurred by the effect of the radiation. 
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TABLE IV: Extended quintessence models: background effects with respect to the corresponding minimally coupled (QCDM) 
models (according to ^ I VI 13 31) . The models are labeled by EQm — f/DNn, the number m defining the inverse power law 
potential, Eq. (I124H . f/DN standing for the (quadratic) non- minimal coupling in Jordan frame or for the exponential Damour- 
Nordvedt type coupling in Einstein frame, and n singling out the specific model, namely the couplings parameters (see text). 
For the models EQ6 — £4 and EQll — flit is zi/ 2 — 0.49 and z^/i — 0.47, respectively; otherwise it is Z1/2 — 0.48. 



model 


Yi lo. 1 

Al/ 2/ Ai/2 


W?, /(W?/o) Q 

1/2/ V 1/2/ 


(D 1/2 /D )/(- ■ ■ ) Q 


Pol Pa 

u / 


-^0 / -^1 / 2 


E06 - f 1 

i ' Vet' yj t _i_ 


1.021 


1.031 


1.006 


1.114 


1.0387 


EQ6 - £2 


1.000 


0.994 


1.000 


1.008 


1.0032 


EQ6 - £3 


1.000 


0.999 


1.000 


1.005 


1.0003 


EQ6 - £4 


1.031 


0.011 


0.989 


0.982 


0.0985 


EQ6 - DN4 


1.005 


1.055 


1.020 


2.101 


1.0081 


EQ6 - DN5 


1.000 


1.002 


1.004 


1.135 


1.0015 


EQll - £1 


1.005 


1.066 


1.044 


1.211 


1.0807 


EQll - £2 


1.000 


0.989 


1.002 


1.006 


1.0060 


EQll -£3 


1.000 


0.999 


1.000 


1.000 


1.0006 


EQll - £4 


1.029 


1.116 


0.980 


0.994 


0.0975 


EQll - DN4 


1.000 


1.003 


1.027 


1.727 


1.0026 


EQll - DN5 


1.000 


1.000 


1.005 


1.099 


1.0004 



TABLE V: Extended quintessence models: maximum deviation on the convergence power spectrum, AP K {£), aperture mass 
variance, A(Af 2 p (#)) , and shear variance, A( , y 2 (9)), from the corresponding minimally coupled models (QCDM). Every weak 
lensing observable is evaluated at two angular scales and the effects of the linear and non-linear regime (values within parenthesis) 
are accounted for. In second and third column, deviation from general relativity parametrized by 7 PPN (z) — 1, Eq. 1171 . is 
quoted. Model labeled as in table EHH. 



model 




( 7 PPN _ 


1) x 10" 3 


AP K (e) (%) 


A(M 2 ap {9)) (%) 


A(7 2 W> (%) 






z = a 


2 = 1 


i = 180 


t = 7200 6 


e = 2° 


9 = 3' 


e = 2° 


9 = 3' 


EQ6- 


£1 


140 


105 


20 


27 (16) 


18 


27 (18) 


12 


19 (21) 


EQ6- 


£2 


1.6 


1.0 


1.8 


2.2 (0.8) 


< 2 


~2(1) 


1 


< 2 (< 1) 


EQ6- 


£3 


1.5 x 10~ 2 





0.7 


0.7 (0.7) 


0.6 


< 0.7 (< 0.7) 


< 0.6 


< 0.7 (0.8) 


EQ6- 


£4 


31 


20 


4 


8(2) 


4 


£ 8 (< 2) 


4(2) 


6(2) 


EQ6- 


DN4 


5.3 


125 


107 


133 (165) 


98 


130 (176) 


79 


110 (154) 


EQ6- 


DN5 


0.27 


0.61 


14 


< 16 (19) 


13 


16 (19) 


11 


14 (19) 


EQ11 


-a 


312 


293 


27 


32 (66) 


18 


63 (35) 


4 


33 (31) 


EQll 


-£2 


4.4 


3.3 


1.8 


3.2 (0.3) 


< 2 


3(< 1) 





< 2 (< 1) 


EQll 


-£3 


4.2 x 10~ 2 





0.1 


0.3 (-0.1) 


0.1 


< 0.2 (-0.1) 


< 0.1 


^ 0.2 (0) 


EQll 


-£4 


67 


51 


4(18) 


11 (15) 


5 (20) 


11 (14) 


(26) 


7(17) 


EQll 


-DN4 


0.32 


0.83 


59 


77 (129) 


54 


75 (133) 


41 


61 (88) 


EQll 


-DN5 


1.6 x 10~ 2 


4.9 x 10~ 2 


9 


11 (17) 


8 


11 (17) 


6 


9 (14) 



"The "old" Solar System bound is | 7 PPN - 1| < 2 x 1(T 3 . The 
"new" Cassini bound is 7 PPN - 1 = (2.1 ± 2.3) X~10~ 5 . 
b £ = 180, 7200 correspond to the angle 9 = 2°, 3' respectively. 



VIII. DISCUSSION AND CONCLUSIONS 

We have investigated the imprint of quintessence mod- 
els, both in general relativity and scalar-tensor theories, 
on weak lensing observations. 

For that purpose, we have derived the lensing quan- 
tities in a way that does not assume general relativity. 
Then, we have developed a numerical extension to our 
CMB code that computes the shear and convergence 
statistics. This allows us to "CMB normalize" all spec- 
tra. 

We have investigated various models and reach differ- 
ent conclusions: 

1. Concerning quintessence models, density perturba- 



tions grow more slowly than in a ACDM model. 
This imply some difference of about 10%-20% in 
the linear regime that can be amplified to more 
than 50% in the non-linear regime, affecting both 
the shape and amplitude of the spectra. This con- 
clusion was reached in various previous investiga- 
tions. 

2. For scalar-tensor gravity we have shown that, given 
the constraints in the Solar System, a non-minimal 
quadratic coupling in Jordan frame will not change 
the prediction by more than 1%. 

3. On the other hand, runaway dilaton models that 
incorporate attraction toward general relativity can 
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FIG. 13: Other possible signatures of the runaway dilaton models, (left) The spectrum of the density perturbation P m differs 
from the power spectrum of the metric perturbation because of the modification of the Poisson equation. The lensing is sensitive 
to the deflecting potential $ and thus enable to measure <5 e ff fEa. I81[l . (middle) The anisotropic stress has two contributions 
(see Eq; IB6II that are important only on large scales, the ansiotropic stress of the radiation and the contribution of the scalar 
field. The latter is dominant by a large factor in scalar-tensor theories, (right) The WMAP data (cross) compared to the 
runaway dilaton model prediction assuming an optical depth of r = 0.16. for m — 11 (dash) and m — 6 (dot) compared to a 
ACDM. 



lead, for the same scalar field potential, to change 
of order 10% in the predictions. Besides the effect 
on the amplitude, there exists a differential effect 
which modifies the shape of the spectra. This opens 
some hope to be able track such a coupling. 

To illustrate the effect on the growth of density pertur- 
bations, one can look at the value of ag, see table lVll Let 
us also note that the redshift dependence of the source 
distribution plays an important role. 

There are however some hypothesis and limitations of 
our analysis that have to be stressed. First, while investi- 
gating the non-linear regime, we have adopted universal 
mappings. These mappings are calibrated on A^-body 
simulations for ACDM models. We have argued that 
these mappings should hold for quintessence models and 
in scalar-tensor theories, as long as we do not enter the 
strong field regime. The verification of this hypothesis 
will require devoted numerical simulation but we do not 
expect the order of magnitude of the effects discussed 
here to change drastically. 

A second point to be mentioned is related to the possi- 
bility that our universe has been reionized. Reionization 
affects the global normalization of the CMB anisotropics 
and thus our normalization of the density parameters. As 
an example of its effect, table I VII compares the results for 
vanishing optical depth (upper) and for an optical depth 
t = 0.16. Generically it changes the value of erg by 10% 
to 20%. But, it will not change the shape of the power 
spectra. The effect of the reionization on the normaliza- 
tion can easily been understood by looking at the CMB 
angular power spectrum (Fig. I14|) . 

Let us also emphasize an effect that may be of impor- 
tance while interpreting lensing data on large scale, as 
could be obtained by a wide field imager. The radiation 
anisotropic stress of the radiation implies, even if gravity 
is described by General Relativity, that the two gravita- 
tional potentials are not equal so that the value of the 
deflecting potential is not equal to twice the gravitational 
potential. 

Weak lensing observables combine the effects of the 



TABLE VI: Root mean square of the variance of the mat- 
ter density contrast on scale of 8/i _1 Mpc today, as, assuming 
vanishing optical depth (upper) or accounting for the reion- 
ization with optical depth r = 0.16 (lower). 



coupling 


m = 


m = 6 


m = 8 


m = 11 


no 


0.90 


0.50 


0.41 


0.32 


ei 


0.95 


0.57 




0.40 


£2 


0.95 


0.50 




0.32 


£3 




0.50 




0.32 


£4 


0.95 


0.48 




0.31 


DN1 


0.83 








DN2 


0.83 








DN3 


0.83 








DN4 




0.71 




0.40 


DN5 




0.53 




0.33 


no 


1.11 


0.58 


0.48 


0.37 


£1 


1.11 


0.67 




0.47 


£2 


1.11 


0.58 




0.37 


£3 




0.58 




0.37 


£4 


1.11 


0.56 




0.36 


DN1 


0.96 








DN2 


0.96 








DN3 


0.96 








DN4 




0.82 




0.46 


DN5 




0.62 




0.39 



background properties and the growth of the perturba- 
tions and extend from the linear to non-linear regime. 
These two regimes are complementary mainly because of 
the sensitivity on the time at which the modes enter the 
non-linear regime. In conclusion weak lensing survey ap- 
pear to be a key observation of the shallow universe to 
investigate the nature of the dark energy. 
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FIG. 14: CMB angular power spectra: the WMAP data 
(cross) compared to two models differing only by their op- 
tical depth t — (solid) and r = 0.16 (dash). This shows the 
influence of the reionization when one chooses to normalize 
to the CMB temperature anisotropics. 
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APPENDIX A: BACKGROUND EQUATIONS 

We summarize the equation of the background for the 
Lagrangian JQ| following the general equations presented 
in Ref. [?|. The metric takes the form 



ds 2 = a 2 (77) [-drj 2 + - , ; d.r'. l.r'^ 



(Al) 



and we denote the derivation with respect to the confor- 
mal time 77 by a dot and we define the comoving Hubble 
parameter by 



The Einstein equations give the two Friedmann equa- 
tions 



3F(H 2 + K) = 8irG*a 2 p+^tp 2 



2F [H 2 - U + K 



+a 2 U - 3HF ilp ip 
8irG*a 2 (p + P) + tp 2 



(A6) 



+ F wV ? 2 + F tV (ip - 2Hip) (A7) 

For completeness, let us give the analog equations in 
Einstein frame where the metric takes the form 



(A8) 



ds 2 = al(i])[-di] + ~/ lj dx l dx J ] . 
The Klein-Gordon equation takes the form 

(p\ + m*tp* = -alV v - 47rG*a(^)(p st - 3P„) (A9) 
and the Friedmann equations become 



3 (H 2 



K) = 8nG,alp. 



+ <£ + 2Va 2 , (A10) 



3H* = -4irG*a 2 (p*+3P* 
-2ipl + 2Val. 



(AH) 



The scale factors of the Jordan and Einstein metrics are 
related by 

a = A{ip«)a* (A12) 
so that the two cosmic times are related by 

d.t = A{y*)&U, dr j = dr l . Jf . (A13) 
This implies that redshifts in both frames are related by 



Ji = a/a. 



(A2) 



Let us start by the conservation equations. The Klein- 
Gordon equation for the scalar field takes the form 



<p + 2Htp 



-a 2 U^ + 3F V f n 2 + n + K 



(A3) 



while the conservation equation of the matter fields is 
given by 



p + m{p + P) = 0. 



We have used the notation 



F 

' v ~ dtp' 



(A4) 



(A5) 



Note that the matter energy density scales as p oc a~ 3 . 
This is one of the reasons for which it easier to solve the 
system in Jordan frame since in Einstein frame /?* = A A p 
will behave as A^a^Jet" 3 which entangles the evolution 
of the density with the one of the background. 



1 



^) 



(A14) 



and that the physical length associated with to a comov- 
ing length £ c are connected by 



^phys _ £phys 



(A15) 



It will be convenient to decompose the energy-density 
of the scalar field as 



Pip = PMC + PF 



with 



Pf a 



\^ 2 +a 2 U, 

-mF v y. 



(Aie) 

(A17) 
(A18) 



Then, we can define the density parameters as 



n f = 



inGcttPfa 2 
W 2 



p v a 
3FH 2 ' 



K 



K 



(A19) 
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so that the the Friedman equations takes the form 

^2n f + ii v + si K = i. (A20) 
/ 

Note that we could have defined 

= _ 8TrG cS ,oPxa 2 ~ _ P v a 2 ^ _ KF 

X ~ 3W 2 ' V ~3F W 2 ' WF - 

(A21) 

The two sets of density parameters agrees today and are 
related by 



n = qf/Fq. 



(A22) 



We deduce that in a matter universe the Friedmann equa- 
tion can be rewritten as 

% = n m ,o$(l + z) + ^,o^-^(l + z)- 2 - Q K , 

(A23) 

To finish, let us define the function 

E(z) - (1 + z)H(z)/H Q (A24) 
that is explicitly given by 



Fn 



F, 



E A = n mfi ^(l + zr + n^4 — - n K( l + zy Z - (A25) 



F 



APPENDIX B: PERTURBATION EQUATION 

The general gauge invariant perturbation equations in 
scalar-tensor that are being integrated in the CMB code 
have been presented in Ref. [9j. We just summarize the 
scalar mode equations in the case of a single fluid in a 
universe with Euclidean spatial sections. 

In Newtonian gauge, the metric takes the form 

ds 2 = a 2 { V ) [-(1 + 20)drr 2 + (1 - 2i/j)dx 2 ] . (Bl) 
The fluid velocity perturbation is decomposed as 



Su^ = a(-<j),d k V). 



(B2) 



where w = Pj p. 

The fluid conservation equation is given by 



5 N Y • w 

= -AV + 3ijj-3H- r (B4) 

1 + w J 1 + w y ' 



while the Euler equation takes the form 



V + HV 



-o 



r + -A* 



(B5) 



1 + w 1 + w 
where c 2 = dP/dp. T is the entropy perturbation and 7f 
the anisotropic stress. 

Among the four independent Einstein equations, we 
can retain 



ip - <j> = 87rG*P7f + -£5(p, 
b 



(B6) 



2F{i]; + H(f>) + F(f> = -8TrG«p{l + w)a 2 V 

+tpStp + SF- HSF, (B7) 



2FA</> + {ip 2 - ZHF)(f> - 3Ftj> 



-F m 



a + 3 (n 2 + ^ 



Sep 



F ■ 



+ (U v a 2 + 3Htp)5<p + tpS(p + 3-5F. (B8) 

F 

The Klein-Gordon equation for the evolution of the the 
scalar field is then given by 



8(p + 2HS<p - 



A + StHFpsp — U v(p a 



Sep 



3ip)5ip - 2a 2 <t>U v 



2iP)+3(i> + 3Hip+H(p) F V .(B9) 



5 is the density perturbation in Newtonian gauge and In all these equations, we have 



we introduce the density perturbation in comoving gauge 
as 



F = F^Sip, F = F w <p 2 + FtpSip, 6F = F^Sip. 



3H{l + w)V 



(B3) 
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